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Abstract 

We present a random matrix interpretation of the distribution functions which 
have appeared in the study of the one-dimensional polynuclear growth (PNG) model 
with external sources. It is shown that the distribution, GOE^, which is defined as the 
square of the GOE Tracy- Widom distribution, can be obtained as the scaled largest 
eigenvalue distribution of a special case of a random matrix model with a determin- 
istic source, which have been studied in a different context previously. Compared to 
the original interpretation of the GOE^ as "the square of GOE" , ours has an advan- 
tage that it can also describe the transition from the GUE Tracy- Widom distribution 
to the G0E2. 

We further demonstrate that our random matrix interpretation can be obtained 
naturally by noting the similarity of the topology between a certain non-colliding 
Brownian motion model and the multi-layer PNG model with an external source. 
This provides us with a multi-matrix model interpretation of the multi-point height 
distributions of the PNG model with an external source. 
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1 Introduction 



The universality of eigenvalue correlations in random matrix tlieory(RMT) has been aris- 
ing in various fields in physics. The most recent appearance is in the one-dimensional 
Kardar-Parisi-Zhang(KPZ) universality class [1] which is a fundamental universality in 
non-equilibrium statistical mechanics. It has been found that the fluctuation of certain 
physical quantities in the KPZ universality class is equivalent to that of the largest eigen- 
value in RMT. This fact was first pointed out in the polynuclear growth(PNG) model [2] 
and the totally asymmetric simple exclusion process(TASEP) [3] based on the Baik-Deift- 
Johansson theorem in random permutations [4]. 

The PNG model is a stochastic surface growth model and is the most well-understood 
model in the relation with RMT. Based on the results on the symmetrized random permu- 
tations [5-7], the properties of the height fluctuation of the PNG model have been analyzed 
for various boundary conditions. In [2] the height fluctuation of the PNG droplet in an 
infinite space was analyzed. It turned out that the height fiuctuation is described by the 
GUE Tracy- Widom distribution [8]. On the other hand, in case of the growth on half line, 
the height fluctuation at the boundary is described by the GOE and GSE Tracy- Widom 
distribution [9] according to the strength of an external source at the boundary [10]. Fur- 
thermore the equal time correlation of the height fluctuation at distinct points has also 
been analyzed for the PNG model. It has been found that it is equivalent to the process of 
the largest eigenvalue in the Dyson's Brownian motion model between GUEs in the case of 
the growth on inflnite line [11, 12] while in half inflnite case, it corresponds to the process 
in the model which describes the transition between GOE/GSE and GUE [13]. 

In addition to the Tracy- Widom distributions for the fundamental three classes in 
RMT, GOE, GUE and GSE, two distribution functions, which had not been discussed in 
the study of RMT, arose in the description of the height fluctuation of the PNG droplet 
on inflnite line with external sources at boundaries [10, 14]. These are called Fq and 
GOE^. These functions also appear in the study of current fluctuations in TASEP [15,16]. 
There are also distribution functions which describe the transitions between two functions 
mentioned above. For instance, there is a distribution function, depending on a parameter 
cu, which approach the GUE Tracy- Widom distribution as a; — > oo, tends to the Gaussian 
as (x^ — > — oo and becomes GOE^ when u = 0. We have not yet known what kind of 
random matrix ensemble corresponded to Fq. The GOE^, on the other hand, can be 
understood in terms of the language of RMT. It is written as the square of the GOE Tracy- 
Widom distribution and can be interpreted as the distribution function of the larger of the 
largest eigenvalues of two independent GOEs. Its Fredholm determinant representation 
was obtained in [17]. There, it was also shown that the scaling limit of certain eigenvalue 
statistics of two independent GOEs is described by its kernel. 

In this paper, we present another point of view of the GOE^ as the scaled largest 
eigenvalue distribution of a random matrix with a deterministic source. The advantage of 
this ensemble is that it can describe the GUE-GOE^ transition in the PNG model, which 
is induced by the change of the external source while this transition cannot be obtained 
by the measure of "the square of GOE" . 
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From the point of view of the PNG model with external sources, our viewpoint is more 
natural than that of "the square of GOE" . We turn our attention on the topology and the 
structure of the determinantal measure in the multi-layer version of the PNG model with 
an external source. Our random matrix ensemble is obtained by noticing that a certain 
type of the non-colliding Brownian motion has a similar topology and the structure of the 
measure. 

The random matrix obtained in this paper is a special case of the random matrix with 
a deterministic source, which has been studied in [18-26]. In a special case, which is 
different from ours, its eigenvalue statistics has an interesting property that the level spac- 
ing distribution becomes an intermediate one between the Poisson and the Wigner-Dyson 
statistics. This property connects it with the fields in physics such as the quantum chaos, 
the metal-insulator transition in disordered systems and so on. Recently the universality 
of such eigenvalue statistics in this model have been analyzed by the Riemann-Hilbert 
technique [27-30]. 

This paper is arranged as follows. In the next section we summarize the results on 
the height distribution in the PNG model with an external source. Then we introduce a 
random matrix with a deterministic source and compare its largest eigenvalue statistics 
with the height fluctuation of the PNG model. In section El we consider a reason why 
such a random matrix model appears in the context of the PNG model with an external 
source. We notice that the multi-layer version of the PNG model can be regarded as a 
non-colliding many body random walks (vicious walk) and hence inherits a free fermionic 
picture. By considering a certain vicious walk which is expected to have a similar property 
to the multi-layer PNG model, we see there appears naturally a Dyson's Brownian motion 
model (a multi-matrix model) involving an external source in the initial condition. Finally 
we show the equal-time correlation of the PNG model with an external source is equivalent 
to the dynamical correlations of the largest eigenvalue in the Dyson's Brownian motion's 
model. The last section is devoted to the conclusion. 

2 PNG model and random matrix with deterministic 
source 

2.1 PNG model with external source 

First of all we give the definition of the one- dimensional discrete PNG model [12]. Let 
r G {■ • • , — 1, 0, 1, ■ ■ ■ } and t G {1, 2, ■ ■ • } be the space and time coordinate respectively. 
The model is a stochastic surface growth model and consists of the rules (a)~(c). Fig.l 
illustrates these rules. 

(a) At time t, a nuclear is generated randomly at r = — t -|- 1, — t -|- 3, ■ ■ ■ , t — 3, t — 1. The 
height A; (= 0, 1, 2, ■ ■ ■ ) of each nuclear is independent of r and obeys the geometric 
distribution with some parameter which is fixed shortly. 
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This is called nucleation. While a nucleation is stochastic, a growth after the nucleation 
given by the rules (b) and (c) is deterministic . 

(b) Each nuclear grows laterally with one step toward right and left at each time step 
and then forms a step. 

(c) When a step crashes against another step, the height of the crashing point is that of 
the higher one. 

These rules (a)~ (c) can be formulated by a single equation, 

h{r, t + l)= max [h{r - 1, t), h{r, t), h{r + 1, t)] + cj(r, t + l). (2.1) 

Here h{r,t) is the height of the surface at time t and at position r, and uj{r,t) means a 
height of a nucleation which is a random variable. 

In this paper we set the random variable ui{r,t) as follows. 

• When t — |r| > and t — r is odd, 

¥[uj{r,t) = k] = 



Otherwise 



- a^g) (a^g) r = -t + 1, 

- q)q^ otherwise, 

cu(r,t) = 0. (2.2) 



Because of the condition that the nucleation occurs only when t — \r\ > 0, the surface 
grows into a droplet shape. The condition that t — r is odd is imposed only for a technical 
reason and is irrelevant to the physical properties. Basically uj{r,t) is a geometric random 
variable with a parameter q. Only on the left edge, the parameter is a.^{^ < « < ;^)- 
Thus the parameter a adjusts the boundary condition as an external source. This situation 
representing ()2.2|1 is depicted in Fig.2. 

In the setting of ()2.1|) and ()2.2|) . we focus on the distribution of the scaled height at 
the origin defined as 

% = 0,t = 2iV)-aiV 
ti]^ = —i , [2.6) 

where a = yz^, d = The limiting distribution has been obtained in [2, 10, 14]. 

It depends dramatically on the value of a. The result is 

fF2{s), for ^<a<l, 
Fi(s)2, for a = 1, (2.4) 



0, for 1 < a < 
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where -^2(5) and -Fi(s) are the GUE and GOE Tracy- Widom distribution respectively [8,9]. 
^2(5) is defined as follows. Let Ai be the largest eigenvalue in x GUE random matrix 
with the measure, 

^-^rNh^dM2, (2.5) 

where M2 is an N ^ N Hermitian matrix. Taking the edge scaling for Ai, 

Ai = v^+^^, (2.6) 
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we define -^2(5) as the limiting distribution of Xi, 

^2(5) = lim P[Xi < s\. (2.7) 

Fi{s) is also defined in the same way as above for N x N GOE random matrix with the 
measure, 

trMi 2 

e — —dMi, (2.8) 

where Mi is a.n N x N real symmetric matrix. Hence the distribution Fi(s)^, can be inter- 
preted as the distribution of the largest eigenvalue in the superimposition of eigenvalues of 
two independent GOE random matrices. In this sense the distribution defined by -Fi(s)^ 
is denoted as GOE^. 

Both F2{s) and -Fi(s)^ can be represented as the Fredholm determinant, 

F2(s) =det[l + /C2(7], Fi{sf = det [I + ]Cug], (2.9) 

where g{x) = —X{s,oo){x)- The definition of the Fredholm determinant with kernel }C{x,y) 
is 

det[l + Kg] = ^ — ■ j dxi- ■ ■ dxkg{xi) ■ ■ ■ g{xk) dei[K{xi,Xj)]\^j^^. (2.10) 

J —00 J — 00 

The kernels ]C2{x,y) and }Ci2{x,y) have the following form, 

POO 

^2(3;,?/)= / da Ai{x + a) Ai{y + a), 
Jo 

POO 

lCi2{x,y)=1C2{x,y)+ki{x) daM{y-a). (2.11) 



Here ]C2{x,y) is called the Airy kernel [31] and fCi2{x,y) is obtained in [17]. 
Thus when a < 1, the result ()2.4|) can be rewritten as 



lim ¥[Hn <s] = det[l + K,g], (2.12) 



where 



^(^^ ^ ) ^2(3;, y) = /o°° daAi{x + a)Ai{y + a), for < a < 1, 

' ICi2{x, y) = IC2{x, y) + Ai{x) da Ai{y — a), for a = 1. 



When 1 < a < -i= , the result ()2.4|) that the hmiting distribution vanishes imphes that 
we have to define another scaled height, 

= hir = a,t = 2N)-aaN^ ^^^^^ 

dcN2 

where ac = '?^(^~^"'?^+" ) an^l do = "^g^^^'^^)^*^" _ Then the limiting distribution 

(a-q2)(l-ag^) (l-agS ) 

becomes the error function, 

lim F[H'ff'^ <s] = ^= \ die-^-, for 1< a < — . (2.15) 



27T J-oo 

2.2 Random matrix with deterministic source 

When we look at the above results, we notice a fact that the distribution at a = 1, which 
we call GOE^, appears in the intertwining point between Gaussian statistics and RMT 
statistics. First, let us consider the distribution of the height at the origin when a is 
extremely large. In this case we can easily expect that it is the Gaussian since the height 
at the origin is virtually determined by nucleations at the only one point, the left edge due 
to the rule (b). Equation ()2.15|) confirms this discussion. When a = 0, on the other hand, 
it is described by the GUE Tracy- Widom distribution as discussed above. The GOE^ arises 
in the situation where these two statistics compete i.e. when a = 1. 

Considering this property, we introduce a simple random matrix model which is ex- 
pected to describe the GUE random matrix-Gaussian transition of the largest eigenvalue. 
We consider a random matrix Hq defined as 

Ho = H + V, (2.16) 

where H is the GUE random matrix and V is the deterministic diagonal matrices, V = 
diag(ej) (j = 1, ■ ■ ■ , A^). The measure can be defined as 

e-'^'^'dHo. (2.17) 

Focusing on the case where 

ei=e,e2 = ■■■ = eN = 0, (2.18) 

we can expect that when e is sufficiently large the largest eigenvalue is isolated from others 
by the effect of e and then the fiuctuation becomes the Gaussian while this becomes the 
GUE Tracy- Widom distribution when e is sufficiently small. 

The eigenvalue statistics of ()2.16p and ()2.17p is analyzed in [18-26]. For the distribution 
of the largest eigenvalue Ai, their results lead to the Fredholm determinant expression, 

P[Ai <s]= det[l + KNg]. (2.19) 
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Here the kernel KN{x,y) can be expressed as a double integral, 



KN{x,y) 



N 



dv f du T-r 

r27ri J 27r-l-l 



'^j 2 



e 4 



-V +iux+2vy 



(2.20) 



where F denotes a contour enclosing {ej}j=i ...^tv anticlockwise and the integration wrt u 
is taken from — oo to oo and not to cross F. 

In this formula with the condition ()2.18|) . we consider the following scaling limit, 



X 



'2N + 



X 



V2Ne 



y 



'2N + 



Y 



A 



(2.21) 



with A(< 1) fixed. This corresponds to the edge scaling limit (j2.6j) of the eigenvalues. 
Applying the saddle point method, the limiting kernel can be calculated. The result 
depends on the value of A, 



/C2(X, Y) = /°° da Ai{X + (T)Ai(r + a) 



for A < 1, 



^ ^/Ci2(X,F) = /C2(X,F)+Ai(X)/;"c/aAi(y-a), for A = 1. 

Note that the kernel ()2.22j) is exactly the same as that in the PNG model with an external 
source (j2.13|) and A, the parameter of the deterministic source in RMT corresponds to 
a, the parameter of the external source in PNG model. Hence we find that the height 
fluctuation in the PNG model with an external source shares the limiting distribution 
with the largest eigenvalue fluctuation of the random matrix with deterministic source 
deflned in ()2.16|) - ()2.18|) . The same distributions appear in the distribution of the largest 
eigenvalue of the non-null complex sample covariance matrices [32]. 
When A > 1 in (|2.21|) . we change the edge scaling (j2.21|) to 



X = AgN2 + BgX, y = AgN2 + BgY, 
where Ag = 775 ~'~ x) = y ^^^r-- The kernel is asymptotically 



(2.23) 



KM{x,y) ^ 1Cg{X,Y) 



r2Nh{y-x) 



for A > 1. 



(2.24) 



G 



Thus we get 



lim P 



Xi-Ag^/N 



B, 



< s 



G 



'2tx 



die 



(2.25) 



This result corresponds to ()2.15p . Note that the factor t'^^'^y ^) in ()2.24|) does not 
contribute to the Fredholm determinant. 
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At last we give a proof of ()2.24j) . The main result in this section ()2.22|) will be derived 
in section 13.41 under more general situation. 

Rescaling the variables in ()2.2n|l under the condition ()2.18|1 such that, 

[N iu [N . /iV , , 

one finds 

K^(x,y) = j-—r / dz / dw —= ^e^^^«('")-^^(^)^ 

^ {2my Jr' J -too w z-Aw- z {V2Nx-V2Agn)w 



(2.27) 



where the contour F' encloses the origin and A anticlockwise and 



2 

fciw) = Y~ ^^gw + ln«;. (2.28) 

K^{x^y) can be divided into two parts, 

K^{x,y)=K%\x,y)^K%\x,y). (2.29) 

Here and correspond to a contour integral of z around A and the origin respec- 
tively and have the forms, 

) ^ V2yVA£^ e-^foW / ^ ^4 ^ (2.30) 

^""^y^- {2myf'^'J_,J'^wz-Aw-z^{V2N.-V-2AaN)^'' ' ^^-^^^ 

where ^ means the contour enclosing the origin anticlockwise. We analyze the asymptotic 
form of each K^^\ (i = 1, 2) by the saddle point method under the scaling (j2.23j) . 

At first we consider K^j^\x, y). We deform the path of integration of w to the one which 
passes a critical point = A of fciw). We set 

w = A+ (2.32) 
and substitute this into the integrand in ()2.30|) to get 

NfcMr^NfciA)-^, 

^(V2Nx-V2AgN)w ^ ^(V2Nx-V2AGN)A^iwiX 33) 
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Thus one finds 



^(1) _ gvWA(y-.)_^^_ (2.34) 

For Kj^ {x,y), we deform tlie patli of the variable z in a way that it crosses another 
critical point = of /g(^) while as for w, we again use ()2.32|) . Hence z is transformed 
to 

1 iW2 

z=-+ ^ , 2.35 



and then we get 



Nfciz) ~ -Nfc ( i ) - ^ 



2 
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,[V2Ny-V2AGN)z ^ ^{V2Ny-V2AGN) ^iw2Y 36) 

Due to ^n^i and (jT^ . one finds 

iriJ)(x,i/)~e^{^^(^)-^^(^)}e^(^^) 

:(A2-J^-41nA)^0(iV^) ^ (2.37) 



_JV|'a2_ 1 



where we use A > 1. 

Thus from (jT^ . (^3^1 and (|23Z|), we finally get ^M>- 



3 Free fermionic picture 
3.1 Multi- layer PNG 

In the above section we discussed the limiting distribution of the scaled height at the origin 
in the PNG model and its dependence on the external source. In the next step, we would 
like to consider the multi-point equal time joint distributions. To do this, we introduce the 
multi-layer PNG model from the time evolution of the PNG model [11,12]. 

The rule of the multi-layer PNG model is as follows. When two steps collide, the lower 
step is absorbed by the higher one due to the rule (c). We recover this absorbed height as 
the nucleation in the layer below. This situation and a typical example of the multi-layer 
PNG model are illustrated in Fig. 3. Note that the first layer of the multi-layer PNG model 
represents the shape of the PNG droplet while other layers record the time evolution of 
the growth. 

If we treat the i-th layer of the multi-layer PNG model as the z-th random walker, the 
multi-layer PNG model can be regarded as a non-colliding random walks called the vicious 
walk. Thus the multi-point equal time correlation of the PNG model, which is our target, 
corresponds to the dynamical correlation of the first walker in the non-colliding random 
walks. Note that we treat the space axis in the PNG model as the time axis in the point 
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of view of the vicious walk. The advantage of the mapping is that as we will see below, 
the measure can be described in the form of products of determinants. This determinantal 
structure is associated with the structure of wavefunction in the A^-body free fermions, i.e. 
the Slater determinant. 

In [33], we could obtain the multi-point equal time joint distributions in the PNG 
model having external sources at both edges by introducing the multi-layer PNG model. 
The result specialized to the case considered here is the following. 

We define the scaled height Hj^ij) near the origin as 

h{2cNKt = 2N)-aN 
Hn[t) = —I + r, (3.1) 

dIS 3 

2 

where a and d are already given below ()2.3|) and c = lii^^lH. We set the parameter of the 
external source as 

a = l-^^. (3.2) 

dNl 

The equal time multi-point distribution function near the origin is described by the Fred- 
holm determinant, 

lim P[if7v(ri) < Si, ■ ■ ■ , Hf^{Tm) < s J 

A''— ►oo 

= det [1 + ICQ] , 

CO ^ m /»oo POO 
fc=0 ni=l>^-°° nfe=l>^-°° 

(3.3) 

where Q{Tj,^) = -X(sj,oo)(0 (j = 1,2, ■ ■ ■ ,m). The kernel is 

POO 

/C(ri,6;r2,6)=/Cr(Ti,ei;r2,6)+Ai(6) / ciAe-(-+-^)^Ai(6-A), (3.4) 

where 

r-e.t. , . ^ . ^ - Jr rfAe-^(---)Ai(^i + A)Ai(6 + A), n > r„ 

y~ Loo "-^^ ''^Al(^i + A)Al(^2 + A), Ti < T2. 

In the case where all = 0, another representation of this distribution has been obtained 
via Riemann-Hilbert method by Baik and Rains [14]. 

The kernel }C2'^^^ is called the extended Airy kernel [34] and the process described by 
this Fredholm determinant is called the Airy process [11, 12]. It describes the process of 
the largest eigenvalue in the Dyson's Brownian motion model in Hermitian matrices. It 
also describes the correlation function of the PNG model without an external source. 

Looking at the correlation function ()3.3p at one point case (m = 1), we can find that 
it expresses the GOE^-GUE transition induced by both parameters of the external source 
u and a position r. Setting = ri = in ()3.4|1 in the one point case, we get /C = K-u 
which is the kernel of the edge scaling in GOE^ ensemble while in the case where u; ^ oo 
or r oo one finds /C ~ /C2 and the GUE case is recovered. 
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3.2 Non-colliding Brownian motion and multi matrix model 



In subsection 12.21 the random matrix with a deterministic source, which corresponds to 
the PNG model with an external source was introduced from an intuitive argument. Here 
we show that the measure of the random matrix with a deterministic source is obtained 
naturally by considering the measure of the multilayer version of the PNG model with an 
external source. The analysis not only reproduces the result in subsection 12.21 but also 
provides a multi- matrix model which has the kernel ()3.4j) in the edge scaling limit. 

As was shown in the above subsection, it is found through the multi-layer PNG model 
that the structure of the vicious walk (free-fermionic structure) is hidden in the measure of 
the PNG model. Particularly we focus on the two properties in the measure of the vicious 
walk. The first property is the determinantal structure which the vicious walk and RMT 
share. The second one is the topology of the vicious walk. In our case, Fig.(3-b), the 
top layer, which corresponds to the top walker in the point of view of the vicious walk, 
seem to start from the point far from the other walkers due to the effect of the external 
source. In the study of vicious walks, it has been known that in some cases the measure 
can be represented as that of a multi-matrix model in an appropriate scaling limit and its 
topology determines the universality classes of the measure in RMT [35,36]. 

Keeping in mind the above two properties, we begin our discussion with writing down 
the measure of non-colliding Brownian motion which can be regarded as the continuous 
limit of the vicious walk. Let x"J\r = 1, ■ ■ ■ ,N,j = 0, ■ ■ ■ ,M+1) be the position of r-th 
walker from top at the time labeled by j. By the Karlin-McGregor theorem, the measure 
that the walkers have a configuration {xi''^ } is obtained as a product of the determinants 
of propagaotrs of the one-body Brownian motion, 

-^U^et[^,,^M'\x(r'^)]l^,. (3.6) 



Z 



Here Z is the normalization constant and iljjj+i{xi''\xi^~^^^) is the propagator of the one 
body Brownian motion from xi''^ to xi'''^^^ with interval Tj 
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^3.^i^\^^'^) = ^^exp{-(xi^-+^) -4^-))V2T,} . (3.7) 
Then we take the following scaling limit, 



J=0, 



Tj = r,T (j=0,l,2,--- ,M), x^^^ = { j = l,2,---,M, (3.8) 

6„ j = M+l. 

Note that if we take the parameter of the initial condition to be ei = e, 62 = ■ ■ ■ = cat = 0, 
the topology of the configuration becomes similar to that in the multi-layer PNG model 
with an external source. (See Fig. 4.) Substituting ()3.8|1 to fl3.7|) . we get 

det [^,,+i(x(^-), xi^-+^))] = ^t[X^'\-^-A^''^] ^ ^3_g^ 
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where Z' is the normahzation constant. The weight Wt has the following form, 

= lllje"^^^^'- > II det[e'=^^'- i ^ 

i=i ''=1 i=i 

where A*^*-* means < Ar*^ 



Afel 



■det 



(3.10) 



\r r and 

=!,■■■ ,7V 



2 V^j-i 



TV 



2ro 



A{x) = n 



(3.11) 



In ()3.10|1 . we include the term A(e) in order to get a finite expression when we consider 
the special case where ()2.18j) in the later discussion. Applying to the last term in Wt the 
following formula, 



lim 



det 


7 


3 


det 







n 



l<i<j<Af 



6 -^j+j-i 
3 -i 



(3.12) 



where {^j} G C^, one finds for large T, 



det 



det 



det 



AM) 



N-j 



n 



e 

/i^(6, - 6,) A(A(^'^)) 



■det 



AM) \ N-i 



e ^ 



l<i<j<N 



N(N-l) • 



(3.13) 



Thus we get 



lim 



4nn 



M-1 

JJ det[e'=^'^'''^^ 



o..j detyo--^-- ^ ^^^^^^^^^ ^^^^^^ 



j=lr=l j=l ^ ^ 

Up to a trivial constant, ()3.14|1 is equivalent to the weight of the eigenvalues of the Her- 
mitian multimatrix model, 



M M~l 



M- 



(3.15) 



Here Hj{j = 1, ■ ■ ■ , M) is an x hermitian matrix and V = diag(e). This equivalence 
can be shown through the Harish-Chandra,Itzykson-Zuber integral [37,38], 



(3.16) 
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where U is NxN unitary matrix, dU is the Haar measure and A and B are NxN Hermitian 
matrices with eigenvalues and bi respectively. Note that in the case of one matrix model, 
(I3.15|l is the same form as the random matrix with deterministic source ()2.16|1 . ()2.17|1 
discussed in section |21 We further transform the variables Sj and rj as 

\ — g2(ij-tj+i) 

Sj = e^'\ rj = — — — , ro = l/2, r^/ ^ oo, (3.17) 

with ti = 0. Then one finds ()3.15p can be represented as the form of the Dyson's Brownian 
motion model, 

exp[-tiHf + iiV Hi]dHi ■ ■ ■ dHu- (3.18) 



M-l 

n 



.7 = 1 



3.3 Dynamical correlation function 

In the above discussion we construct the Dyson's Brownian motion model with the de- 
terministic source from the non-colliding Brownian motion with the topology appropriate 
to the multi-layer PNG model with an external source. Next we would like to analyze 
the process of the largest eigenvalue of this model since it corresponds to the equal time 
correlation of the height distribution in the PNG model. 
The measure for the eigenvalues in ()3.18|) is 



M-l 



ndet[0(t,,A(^-);t,+i,Ai^-+^))]^^,=ie 



(1)^ 



det 



N 



det 



j,k=i 



Af") 



j=l,...,N 
k=0,...,N-l 



where we denote the eigenvalue for the matrix Hj as 



(3.19) 



4>{U,x;tj,y) 



7r(l— e' 



2(*i-t,) 



exp 



for ti < t 



for ti > tj. 



Under this measure, we focus on the following probability. 



(3.20) 



P 



(A/) 



(3.21) 



We analyze this quantity by following the strategy of Johansson [12] although we can 
calculate it also by the method of Brezin-Hikami [19] with some effort. 

In [12] it is shown that the correlation function ()3.2ip can be expressed as the Fredholm 
determinant for the case where the measure is given in the determinantal form such as ()3.19|) 
by generalizing the method of Tracy- Widom [39]. Applying the procedure to ()3.19p . we 
get 



P 



det [1 + Kg] 



(3.22) 
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where g(tj,x) = —X{sj,oo){x). The kernel is 

K{tr, x; ts, y) = K'{tr, x; t,, y) - (piU, x; ts, y), (3.23) 

where 



K'{tr,x;ts,y) = 


N-1 

^ '^j{x,tr;tM){A~^)j^k^k{to;'ts,y), 

j,k=0 


(3.24) 


^!j{x,tr]tM) = 


POD 

/ 4){tr,x;tM,y)y^dy, 

J — oo 


(3.25) 


<^j{0;ts,y) = 


POO 

J — oo 


(3.26) 


Aj,k = 


POO 

/ dxdye''+'''e-''^(l){0,x;tM,y)y'', 

J —oo 


(3.27) 



The analysis of the kernel using above equations seems to be difficult since ()3.24|1 involves 
the inverse matrix. However we can tackle this difficulty by improving the method of 
orthogonal polynomial in RMT [40]. We use efficiently the degree of freedom that the 
value of a determinant is unchanged by elementary transformations. Let us define Fk^^{x) 
by 

F^x) = k\2'/^ / -, (3.28) 



[z - e 



where r(e') represents the contours enclosing all points (Z = 1, ■ ■ ■ , + 1) anticlockwise. 
We also define Gk,e{x) by 

nk/2 p ^ 

= / t^we'^'/^-VI-- _ e;)^ (3.29) 

where = ei/\/2 (/ = 1, ■ ■ ■ ,N) and 7 represents an arbitrary path running from —ioo 
to ioo. The representations of Fj^^[x) and Gj,e(a;) in ()3.28j) and ()3.29|1 corresponds to 
the multiple Hermite polynomial of type I and II respectively discussed in [28] except 
some prefactors. Note that when = 0, both ()3.28|) and ()3.29|) become the integral 
representation of the Hermite polynomial with degree k. 

One can easily check that Ffc(x)'s (resp. G^'s) are linear combinations of e'^'^^'s (resp. 
x'^'s). Hence one finds 

det[e^^^^'^]ffc=i = const, x det[F,- ,(Ai'^)]j=o,...,jv-i, (3.30) 

k=l,...,N 

dei[{\f''^f] j=i,...,N = const. xdet[Gke-^Me{^f^)] j=i,...,n . (3.31) 

fc=0,...,Af-l ' k=0,...,N-l 
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Applying the same procedure as above, one can show that the operator K on the right 
hand side in ()3.21|) can be replaced by another operator K, with the kernel K(tr, x; tg, y), 



K{tr, x; ts, y) = K\tr, x; t,, y) - (piU, x; ts, y), (3.32) 



where 



N-l 



K'{U,x;t„y) = ^j{x,tr;tM){A-^)j,k^k{0;ts,y), (3.33) 

j,k=0 

/oo 

F,-,(x)e-^'V(0,x;t„y)rfx, (3.34) 
■oo 

/oo 

(j){tr, x; tM, y)Gj^^-t,t^{y)dy, (3.35) 
■oo 

_ 2 

Aj,k= dxdyFj^,{x)e (f){0,x]tM,y)Gi,^^~tM^{y). (3.36) 



Substituting (|3.2(jp . (|3.28|) and (|3.29|) into these equations, we get 



^,{0;tg,y) = e-(^+i/2)*^e"^V,.e-...(2/), (3.37) 
^,{x, U; tu) = e'^^+'/^^^''''-''-^G,,,-tr,{x) (3.38) 



and 



/oo 
dx F,,,(x)Gfc,,(x)e-"' = e-(^+^/2)tMy^2^-^-!5,-fc. (3.39) 
-oo 

We chose F and G in such a way that Aj^k becomes diagonal and can hence be easily 
invertible. Note that when = and tM = 0, ()3.39|1 represents the orthogonality 
of the Hermite polynomials. Our treatment here generalizes the method of orthogonal 
polynomials for the multimatrix model between GUE [41,42] to the case with an external 
source. 

Substituting (UniD- (USHD to (UnSl), one gets 

K'{tr.,x;ts,y) = ^j(x, t,,; tM)$j(0; t„ y) 

3=0 

N-l 

= e-^^ 5^G,-,-..,(x)F,,,-.,(y)e-(^-+^/^)(*-*^) 

j=0 

(3.40) 
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where in the last equahty we use ()3.28|) and ()3.29|) . The summation term in the integrand 
can be deformed to 

nil e'^w - t\ 



N) 



N 

= TT ( '^ '^-^l \ ) _^ c-^-tj^ ^ ^ yc 

ii \^ ^ts-, _ y I gtr^ _ g'^ (e^^w — e'jY_ J(e*'-w — e'jy) {e^^w — e'l) ■ ■ ■ (e*'-w — e'^) 

"/e^W^ ' 1. (3,41) 

±J. \^ gis^ ~ ^'i J 1 — e*^2; -^-^ \e*''w — e^/ e*''w — e^^'z J 

From ()3.40|) and ()3.41|) one finally finds the double integral formula of the kernel, 



-y 



{2TTi) 



[ dz [dwT] . \ g-V2-V2^x.-.V2+V2..^ 



(3.42) 



where we use the fact that the second term in ()3.4H) does not contribute to the double 
integral in ()3.40|) . In the one matrix case where tr = tg, the kernel in ()3.42p becomes 
equivalent to the kernel (|2.2Uj) obtained by Brezin-Hikami except the prefactor irrelevant 
for the determinant. 

3.4 GOEVGUE transition 

We consider the asymptotic form of the kernel ()3.23|1 under the edge scaling. Hereafter we 
consider the case where ei = e, e-,- = (j = 2, A^). Rescaling the integration variables 
w,z, the kernel ()3.42j) is rewritten as 

X [ dz f ix ' ~ ^ e^^v-^N. 

Jr'^ w e*22 - e/V2N we*i-*2 - z eV2Nwx-2Nw 

(3.43) 

where the contour F'^ encloses the origin and anticlockwise and 

f{w) = — -2w + \nw. (3.44) 
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We analyze the asymptotic behavior of ()3.43|) by applying the saddle point method. We 
scale the variables in ()3.43|) as 

X = V2N + J\ , y = V2N+ 



^ = v^2iV(l--^). (3.45) 



The critical point of the function f{w), 

w, = 1, (3.46) 

turns out to be the double critical point, 

fiw,) = f'iw,) = 0. (3.47) 

Thus by deforming the paths of the integration in ()3.43|) in a way that they cross the 
critical point, 

^ = 1 - ^, (3.48) 
TVs 

. = 1 + ^, (3.49) 

we get 

where we use 

fiw) ~ fiw,) + f'{w,){w - w,) + ^^^{w - w^f + ^^^^{w - ^c)^ 

= — + —wt (3.51) 

2 3A^ ^ ^ ' 

For other terms in ()3.43p we obtain the asymptotic forms under the scaling ()3.45|) . 

e*^w — e/ \/2N Ti + u — iwi 



e*'^z-e/V2N T2 + u + iw2' 
1 N's 



yjQti-t2 _ ^ — Tl + IWl + IW2 



(3.52) 
(3.53) 



e 



/2Nzy-2Nz 



1 



qV2Nwx-2Nw 
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Hence we eventually finds 



dwi /"°° dw2 ( 1 

H ■ — — I e 



, 27r 27r V r2 - n + i[w\ + W2) + r2 + 2^2 



«(?H«l+6w2) + §(M'i+1il2) 



/•OO /'OO 

X / dAe-^(^i-^2)Ai(ei + A)Ai(6 + A) +Ai(ei) / rfAe-('^+^2)^Ai(6 - A), (3.56) 
where we use the integral representation of the Airy function, 

Ai(x) = / rfAe'^"+t^'. (3.57) 



It is also necessary to consider the asymptotics of (p{ti,x] t2, y) in ()3.2U|) under the same 
scaling as in ()3.45|) . One has 

1 ArV3 f 1 X 1 , x2 

(^2-ri) + — -7^(r2-ri) + ... 



I_e2(ti-t2) 2(r2-ri) [ A^Vs^^ 3iV2/3 

11 7 

+ - 6)(r2 - ri) - ^^(6 - 6)(r2 - ri)2 + Y^^(r2 - r^f 

+ i]^fe-ei)^ + J7^(r2-r0^}. (3.58) 
Substituting these into ()3.20|) . one gets 

^ _ r v/2Ari/6e^^/='(n-r2)+7vV3(^,_^,) joo^ rfAe-^(-i--^)Ai(ei + A)Ai(G + A) for n < r^. 
1 for Ti > r2. 

(3.59) 

Thus we finally find 

^N^l^k e^'''(^-^^)+^'''«-«^)^(ri,ei;r2,6) (3-60) 

where /C(ri,,^i; T2,C,2) is the same as the one appearing in the analysis of the PNG model 
with external source ()3.4|1 . Since the prefactor is irrelevant to the value of Fredholm 
determinants, we have shown that the fluctuation of the PNG height described by ()3.3|) - 
fl3.5|) is equivalent to the fluctuation of the largest eigenvalue of the multi-matrix model 
fl3.18|) with an appropriate choice of V. 
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4 Conclusion 



We have presented a random matrix model with a deterministic source, the largest eigen- 
value of which describes the height fluctuations of the PNG model with an external source. 
Depending on the value of the deterministic source, the distribution of the largest eigen- 
vahic in this model becomes the GUE Tracy- Widom distribution and the GOE^. It also 
describes the transition between these two, which could not be understood by the previous 
interpretation of the GOE^ as "the square of GOE" . Our model gives not only another 
representation of the GOE^ but also a unified picture including GUE, GOE^ and their 
transition. 

We have also considered a vicious walk model which has a similar topology to that of 
the multi-layer PNG model. In an appropriate limit, this leads us to a Dyson's Brownian 
motion model with a deterministic source as an initial condition. At last we have found the 
process of the largest eigenvalue in the edge scaling of this model describes the multi-point 
joint distributions of the PNG model with an external source. 
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Figure Captions 



Fig.l: Rules of the discrete PNG model. Here we draw the solid lines using h{[r],t). 

Fig. 2: Typical situation in the PNG droplet with an external source. In this model the 
rate of growth at the left edge(r = —t + 1) is higher than that at other places since the 
parameter of a nucleation at the left edge, a^^q, is larger than that at other points, ^/q. 

Fig. 3: (a) Rule of the multi-layer PNG model. Due to the rule(c) of the PNG model, 
the part illustrated by the thick lines vanishes in the first layer. This part, however, is 
recovered by the nucleation of the layer below. For the second and subsequent layers, we 
follow the same procedure. 

(b) A typical configuration of the PNG droplet with an external source for a = l,q = ^ 
and t = 200. 

Fig. 4: Non-colliding Brownian motion in the scaling limit (|3.8|) . 
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